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We examine a model of classical deterministic computing in which the ground state of the classical 
system is a spatial history of the computation. This model is relevant to quantum dot cellular 
automata as well as to recent universal adiabatic quantum computing constructions. In its most 
primitive form, systems constructed in this model cannot compute in an error free manner when 
working at non-zero temperature. However, by exploiting a mapping between the partition function 
for this model and probabilistic classical circuits we are able to show that it is possible to make this 
model effectively error free. We achieve this by using techniques in fault-tolerant classical computing 
and the result is that the system can compute effectively error free if the temperature is below a 
critical temperature. We further link this model to computational complexity and show that a 
certain problem concerning finite temperature classical spin systems is complete for the complexity 
class Merlin- Arthur. This provides an interesting connection between the physical behavior of certain 
many-body spin systems and computational complexity. 



PACS numbers: 05.20.-y,05.10.-a,05.70.Fh,89.70.Eg 



2 



The creation of massive digital and deterministic computing systems represents one of the greatest triumphs of the 
last century. Increasingly, however, as the components in our information processing devices shrink to atomic scales pQ, 
the two defining characteristics of these systems — that they are digital as opposed to analog, and that they are 
deterministic as opposed to probabilistic — are beginning to be revealed as approximations that arise from considering 
macroscopic or mesoscopic physical systems. As computing elements are made smaller they are increasingly subject 
to noise and an inability to be perfectly controlled [2j- Even further down this route lie quantum computers (3H7] 
where, instead of probabilistic time evolution, the computer is made of elements which function according to the laws 
of quantum information. These trends point toward the necessity of studying computing elements in settings where 
digital and deterministic functionality are not a priori guaranteed [5] . Here we consider a physical method for classical 
computing wherein a computation is encoded into the ground state of a classical many-body system of spins. At zero 
temperature in this model, computation is deterministic and thus we label this model classical ground state spin 
computation. The model has many predecessors, include quantum-dot cellular automata pHll], the broadcast model 
on trees [TSJ , quantum ground state computing [T3J [2] , and universal adiabatic quantum computing [T5HT!?] . We focus 
our work here on this model when the system has thermalized at a non-zero temperature. At non-zero temperature, 
if we do nothing to compensate for the effects of thermal equilibrium, ground state spin computation fails. Here we 
show that it is possible to use the ideas of fault-tolerant classical computing to design ground state spin computers 
which function deterministically with high probability at non-zero temperatures below a critical temperature. We 
achieve this by mapping the thermal ensembles of our computational spin systems onto the distribution of ensembles 
produced during classical probabilistic computations. This is the main result of the paper: some classical spin systems 
in thermal equilibrium can be thought of as enacting classical probabilistic computations spatially across the system. 

Finally, our mapping leads naturally to problems which are complete for the complexity class promise Merlin- 
Arthur [201 HI] , a result which can be regarded as a finite temperature version of the Cook-Levin theorem 123] 
from the theory of computational complexity. The new computational model we consider thus serves as a bridge 
between statistical mechanics and the classical computational complexity of probabilistic computation. We also 
obtain a computational complexity result concerning the difficulty of identifying when a system admits a mapping 
to a probabilistic circuit by showing that a restricted version of this problem is NP-complete. This implies that the 
problem of identifying when a physical system can be seen to be performing a computation is itself computationally 
intractable. 

The outline of our paper is as follows. In Section [I] we introduce the model of computing with spins in the ground 
state using a particular energy function for this system. In Section[H]we show that at non-zero temperature the model 
from the previous section fails to properly compute. We do this for an extremely simple and previously explored model, 
but present the result using two different methods, one related to reinterpreting the transfer matrix and the other 
related to classical circuits applied to simple thermal ensembles. Generalizing these methods allows us to discuss the 
energy functions arising in Section [I] as probabilistic circuits. This more general formulation we discuss using two 
different methods in Section [III| Motivated by the mappings discovered in Section [III| we then return to our original 
model and define the broader class of energy functions consistent with these models in Section [V] At this point we 
also point out how our models differ from classical models of quantum-dot cellular automata. The more general 
setting leads us to ask questions about how algorithmically hard it is to decide if a given physical system supports 
computation in these models. We show that this problem, even in a very weak form, is NP-complete and thus likely 
to be intractable. In Section |VI| we discuss how to design ground state spin models that, unlike the generic case, 
do compute fault-tolerantly at non-zero temperature. Our construction is intimately related to the original fault- 
tolerance construction of von Neumann [8] and we show that von Neumann's threshold for fault-tolerance is related 
to a critical temperature in our system. Finally in Section VII we discuss how the model we consider is related to the 
Cook-Levin theorem from computational complexity and show how this leads to certain natural problems about our 
model being Promise-MA-complete. 



I. CLASSICAL GROUND STATE SPIN COMPUTING 



Here we introduce the model of classical ground state spin computing. This model, in particular restricted cases, 
is relevant to a variety of different models considered elsewhere. For instance, quantum dot cellular automata can be 
partially modeled by classical ground state spin computing 9, 10J. However the model we consider is more general than 
these specific instances. This larger generality comes along with certain unphysical assumptions: our models contain, 
for example, three-body interactions which are not necessarily easily achievable in a physical device. Physically we 
imagine that models such as the one we consider might emerge in certain limits where effective many-body interactions 
can emerge. However, irrespective of the physical implementation, the model provides a set of classical many-body 
interacting systems which can be mapped onto probabilistic classical circuits. This connects statistical mechanics 
with classical probabilistic circuits thus bridging computer science and physics in a new and interesting manner. 
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We further note that a quantum version similar to this model has been studied by Mizel et al O Q3] • This later 
model contains significant difference owing to the quantum nature of the computation: in particular the ground 
state does not contain the computation laid out spatially, but instead exists in a superposition of the computation 
being carried out spatially. These quantum models are also connected to universal adiabatic quantum computing 
schemes [TSHT!?] and piecewise variations on these schemes 24 27 . By studying the classical analogy of these schemes, 
we hope to shed light on how these later quantum methods can be made fault-tolerant. 



A. The Model 



A combinatorial circuit, C = (G,L), is a directed acyclic graph, G = (V,E) with vertices V and edge set E, where 
the internal vertices, V in , of the graph are logical gates, L : V in — > Q (Q is the set of all logical gates), and the external 
vertices, V ex , are the inputs and outputs to the circuit. External vertices that have no edges leading to them are 
inputs and external vertices that have no edges that lead away from them our outputs to the circuit. If a circuit has 
n input vertices and m output vertices, then to such a circuit we can assign a function / : {0, 1}" — > {0, l} m which 
results from propagating the input through the logic gates to the output vertices. Note that at this point we require 
that fan-outs in our circuit are represented by gates and do not allow circuits which have fan-in gates. A £;-fan-out 
gate is the boolean function fk : {0, 1} — > {0, l} fe given by fk(x) = (x,...,x). 

Consider the following physical system. For each edge of a combinatorial circuit associate a single subsystem with 
only two possible states, and 1, i.e. associate a bit to every edge. Consider next a logic gate which has incoming 
edges labeled by bits i\, . . . , ij and outgoing edges labeled by t\, t^, ■ ■ • , ifc. For such a logic gate, I, define the 
following energy function 

• -{i ffi *>-<'..<*-.« , (1) 

where /; is the function the logic gate is computing. Note that this energy contributes when the inputs and outputs 
follow the rules of the logic gate, but is A otherwise. Now define the energy for a particular configuration of the bits 
associated with edges as the sum over all logic gates of these energy terms: 



E c (s) = J2 E l ( 2 ) 



where each Ei acts on the appropriate input and output bits and s is the collection of bits for the entire system. This 
energy is a function of the labels on all of the edges. Its zero energy configurations are ground states which consist of 
configurations which correctly compute the logical function corresponding to the circuit C. Note that at this point 
the ground state is degenerate: every valid input (and corresponding output) defines a valid zero energy configuration. 
We call an energy function constructed in this fashion from a circuit a circuit energy function. 

Suppose further that, in addition to imposing energy constraints to enforce logic gates, we also impose energy 
constraints to fix a particular input to the circuit. Suppose that the input vertices are Vi,t>2, ■ ■ ■ ,v n and we wish to 
force the input x\, X2, ■ ■ ■ , x n , where Xi £ {0, 1}. Then to each edge whose bit is labeled s, lead away from vertex m 
we can associate an energy term E(si) = if Sj = Xi and E(si) = A otherwise. Taking a sum over all of these terms 
for the input vertices (and corresponding edges) we can thus add a term such that the ground state configuration 
consists now only of the input x\, . . . ,x n propagated to the output of the circuit, f(x%, . . . ,x n ). In particular for 
input x £ {0, 1}™ define the 

n 

EcA*) = E C (s) + J2 E *i( s i) = E c(s) + E x {s), (3) 



where 



„ , , ( if Sj = Xi where Sj corresponds to input Vi , ... 

1) \ A if Si ^ Xi where Sj corresponds to input m 

We call an energy function made up of a circuit energy function plus an input forcing energy function a computed 
circuit energy function. 

The above description of how to take a combinatorial circuit, C, plus its input, (x±, . . . ,x n ), and converting it 
into a many-spin energy function is what we term the classical ground state spin computing (CGSSC) model. Just 
as in quantum dot cellular automata [51 IIP) , the computation occurs when the system is in its ground state and is 
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spatially spread out across the device. Unlike in quantum dot cellular automata, however, this model can implement 
logic gates by using interactions beyond just pairwise interacting spins (or pseudo-spins in the case of the quantum 
dot configurations.) For example, constructing a quantum wire in the CGSSC model will correspond directly to an 
identical model in the semi-classical limit of quantum dot cellular automata models. Gates in quantum dot cellular 
automata, however, are implemented in a way which is not directly analogous to the CGSSC model. In section |V| we 
return to this issue and define a more general set of energy functions wherein our results still hold and compare this 
with quantum cellular automata models. 



B. Example 





FIG. 1. Left: An example circuit whose circuit energy function is given in the text. The bits of the ground state spin energy 
function are labeled Sj. Right: the spin system corresponding to this circuit. Here the spins connected by a triangles have 
interactions between them as specified in the main text. 

To be explicit, let's consider an example circuit involving four inputs, a few logical gates, and one output as 
diagramed in Fig. [l] The energy function for the AND gate is then, for example, 

Eand(si,S2,s 3 ,s 4 , 85) = A((I - si)(l - s 2 )s 5 + (1 - s 1 )s 2 s 5 + si(l - s 2 )s 5 + sis 2 (l - 85)) 

= A(s 5 + s x s 2 - 2s i s 2 s 5 ). (5) 

The full circuit energy function is given by 

OR 



Ec(si,S2,S3, S4,, s 5 ) = A(s 3 + S4 + s 6 - S3S4 - 2s 3 s 6 - 2saSq + 2s 3 S4S 6 ) 

+ A(s 5 + sis 2 - 2sis 2 s 5 ) + A(l - s 7 - s 5 s 6 + 2s 5 s 6 s 7 ) 



(6) 



AND NAND 

Suppose that we wish to force the input to be s% — s 2 = s 3 = and S4 = 1. Then we would add the term 

E C ,x =E C + A( Sl +s 2 + s 3 + (1 - s 4 )). (7) 

II. UNPROTECTED GROUND STATE SPIN COMPUTING FAILS AT FINITE TEMPERATURE 



Above we have defined an energy function whose ground state deterministically carries out a circuit. Having defined 
this energy model we can now consider the physically important question of what happens to this model when the 
physical system described by this energy function is in thermal equilibrium at a finite temperature. Thus we are led 
to consider the Boltzmann distribution corresponding to a non-zero temperature version of a circuit energy function. 
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It then makes sense to consider whether the conditional probability of output f (11,12-, ■ 



(h,t 2 , ...,t m ) given 



the forced input (ii,%2, ■ ■ ■ ,in) f° r this circuit is large enough to distinguish this output from a completely random 
outcome. It is easy to see that for at least some circuits C and inputs x, ground state computation will fail to correctly 
evaluate the function corresponding to the circuit as the size of the circuit being implemented grows. This follows 
directly from examining one of the most basic models in statistical physics, the one-dimensional Ising model |28j . 
This was pointed out in the context of quantum-dot cellular automata in [10] (see also [HIES])- Here we reproduce 
this argument presenting it in two different forms. We do this not just to be pedagogical, but also because these two 
methods will generalize from the one dimensional case to more general circuit computing energy functions. 

Consider the circuit computing energy function corresponding to inputing the bit into a series of n identity gates: 



E I (si,. . . ,s n+ i) = A6 S 



(8) 



where represents the negation of s: — 10 = 1 and -il = 0. It is convenient here, and in the sequel, to work with 
{+1,-1} valued variables instead of the {0,1} valued s^'s. Thus we will define uppercase spin variables to be ±1 
valued versions of the Sj's via Si = 1 — 2si. Then the above energy function can be written as 



A A " 

. . . , S n+1 ) = —(1-Si) + — £(1 - SiS i+1 ) 



(9) 



which we recognize as the one dimensional Ising model with ferromagnetic couplings and a boundary term which is a 
local field. The ground state of this is simply all Sj = (Si = +1), i.e. the initial has been copied by identity gates 
down the line. 

The thermal ensemble arising from this energy function is given by 



Pr(S%, . . .,5 n+ i) = 



eM-PEi(S)) 



where Z is the partition function, 



Z= e*p(-0Ei{S)), 
se{+i,-i}"+ 1 



(10) 



(11) 



and P = (fc^T) -1 is the inverse temperature. It is well known that the one dimensional Ising model does not order 
at finite temperature in the thermodynamic limit [28 . From our perspective, we observe that the system will fail to 
correctly transmit the down the line at finite temperature except in a window of size fc^T < A which goes to zero 
as the system size goes to infinity. 



A. Transfer Matrix Reinterpreted as a Probabilistic Circuit 



The standard method for solving this model is the transfer matrix method which we review and then reinterpret. 
To the energy function add a term dependent on a variable 7 for the last spin: 



E' I (S,j) = E I (S)+ 1 S n+1 . 



1 dhiZ' 

1 di 



Define the two by two matrices 

[Mi] Si+1 ,Si = exp 
where b = /3A, as well as the column vector 

[v]s 1 = exp 



7=0 



1 dZ' 

JZ>lfy~ 



7=0 



:(1 — SiSi+i) 



-2(1-*) 



(12) 



If we calculate Z' for this energy function, we can then use this to calculate the probability that S n +\ is +1 via 
Pr(S n+1 - +1) = (|(1 + S n+1 )) = |(1 + (S„+i)) where 



(13) 



(14) 



(15) 



6 



and the row vector, 



[w T ]s„ +1 = exp[-/?7S n+ i] 



Explicitly the partition function Z' is 



z 1 = ex p 

se{+i 

This can then be written as 



i=l 



Z' = l wT ]s n+1 YllMiWuSiVS! = w T M n M n -i . ..M lV , 

Se{+l,-l} n + 1 *=n 



(16) 



(17) 



(18) 



where the last equation is a row vector, a matrix product, and a column vector thus giving a scalar. This is the 
transfer matrix form of the solution: the partition function is written as a product of "transfer matrices" and, in this 
case, sandwiched between an "initial" and "final" vector. 
Define the following matrices and vectors: 



[Pi]s '+- S < " 1+cxpHO 



Ms, 



We can rewrite Z' as 



[Pm]s - ~ l + exp(-&)' 
Z' = (1 - cxp(-&))-( n+1 ) U ; T P n P„_ 1 . . . P lPm . 



(19) 
(20) 

(21) 



It is at this point that we begin to see our reinterpretation of the transfer matrix method. In particular we note that 
Pin is a vector of probabilities (that is it sums to unity) and P, is a stochastic matrix (its columns sum to unity). 
In other words encoded into Z' is a classical preparation of a probabilistic bit, followed by an evolution according to 
probabilistic gates. Suppose that we let p out = P„P„_i . . . P\p% n denote the output of this probabilistic circuit. Then 



Z' = (1 - cxp(-6)r-» ? 



(22) 



where the sample is taken from the output of the probabilistic circuit. Notice that at 7 = 0, to is a vector made up 
of all components equaling 1. In this case, (w) Pout | = 1 since then this term is a sum over all the outputs to the 
probabilistic circuit. Further 



1 dZ' 



= (i-cx P (-&)r+» 



7=0 



where a = (+1, -1) T . Thus 



(23) 



(24) 



In other words the expectation value of the last spin in the chain is equal to the expectation value for the output of 
the circuit starting with p in and then applying the probabilistic gates Pi, . . . , P n . Note that the combinatorial factors 
(1 + exp(— 6)) 11+1 have canceled out. 

The above description tells us that we can think about the one dimensional Ising spin chain with a forced boundary 
term as a probabilistic circuit related to our ground state spin computation. In particular the probabilistic circuit 
starts with the possibility of an incorrectly initialized boundary. This is formulated in the preparation probability 
vector pi n . Then with probability ilcxp(-h) the ground state spin computation of an identity gate is replaced by a 
bit flip gate. The probability of the final spin of the chain being in a particular state is then directly given by the 
probability that the probabilistic circuit outputs a particular value. 

In this model the probability of the ground state properly computing the desired identity circuit is such that the 
information is washed out at finite temperature. To see this we must actually calculate (S n +\). This can be done 
most easily by diagonalizing the gates. In particular if we define the Hadmard matrix, 



H 



1 

71 



1 1 

1 -1 



(25) 
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then 



Pi = HDH~ = HDH 



(26) 



where 



D = 



1 

n l-cxp(-b) 
U l+exp(-6) 



(27) 



Thus 



(S n+1 ) = a T HD n H Pv 



(28) 



which can be reduced to 



(Sn + l) = 



fl-e- b \ n+1 


' , b~ 




tanh - 


\l + e- b ) 


2 



n+l 



(29) 



From this expression we see that the probability that the final spin is aligned with the input quickly drops to 1/2 
as a function of n (unless b > n, i.e. in a small window of temperature below ~.) Thus the circuit fails with high 
probability as we make the chain longer and longer. This is exactly the result of In this respect, ground state 
spin computing in this model is not tolerant to errors arising from working at finite temperature. Note that we do not 
consider scaling of the temperature to stay within the small window where the computation is correctly performed to 
be a fault-tolerant method as it requires unreasonable physical resources as the system size gets larger. 



B. A Second Approach Using Controlled- Not Gates 



Having shown in the previous subsection how one can reinterpret the transfer matrix method as a probabilistic circuit 
for the one dimensional Ising model with forced boundary, we next present a second way to derive this observation. 
This method is a direct variation on the method used in [30] to study Ising models on Cayley trees. 

On the n + l spins, define the controlled-not on positions i and i + 1 as the function, C,, which maps {+1, 
to {+1,-1}" +1 via 

Ci(S\, . . . , Si, Si+i, Si+2 ■ ■ ■ , S n ) — (Si, . . . , Si, Si+iSi, Si+2) • ■ ■ i S n ). (30) 

We call this a controlled-not because it functions as a deterministic gate which, controlled on the spin Si either does 
nothing to the spin in position (i + 1) (if Si = +1), or flips the spin in position (i + 1) (if Si = — 1.) Further define 
the function which is the composition of controlled-not 's starting at the first spin and working forward (note that 
Ci o C,; + i 7^ Ci + i o Cj): 

C = C n oC n - 1 o---oC 2 oC 1 . (31) 

Notice that if we start with spin 1 in the state Si and all of the other spins S2, ■ ■ ■ , SV1+1 in +1, then the action of C is 
to copy Si down the line: C(Si, +1, . . . , +1) = (Si, Si, . . . , Si). In this sense C is the operation of applying identity 
gates as information is propagated down the spin chain. Note that the Ci are bijections as is C and that Ci is its own 
inverse. Thus 



C- 1 =CioC 2 o...oC„_ioC„. (32) 



Explicitly, 



n+l 



C(Si, S2, S3, . . . , S n+ i) — (Si, S1S2, S1S2S3, . . . , Y[ Si) (33) 



i=l 



and 



C (Si, S2, S3, . . . , S n +i) — (Si, S1S2, S2S3, . . . , S n S n +i). (34) 



Consider the energy function onn + 1 spins: 

£'(S^,...,S; +1 ) = !^(1-S0. (35) 

1=1 

Then clearly a system in the thermal state described by this energy function has each spin in +1 with probability 
1—p = l/(l + exp(— b)) and —1 with probability p = cxp(— 6)/(l +cxp(— £>)). Suppose that we apply the C function to 
this system. That is consider the above system in its thermal ensemble and consider physically applying the controllcd- 
not gates. Then certainly this can be thought of as a system in which the first bit is prepared in with probability 
1 — p, with probability 1—p the controled-nots successfully copy the bit down the line, and with probability p the 
output of the controlled-not is flipped. This is exactly the probabilistic circuit as described in the previous section. 

How is this related to our original El Let the unprimed variables be the spins after C has been applied to the 
system. Then 

Pr(S u ...,S n+1 ) = Pr(S 1 ,...,S n+1 \S[,...,S / n )Pr(S[,...,S' n+1 ) (36) 

s;,...,s; +1 e{±i} 

where 

Pr(5i, . . .,S n+ i\S[, ...,S' n ) = <5c(s;,...,s; +1 ),(Si,...,s„ + i) = ^s{,...,s'„ +1 ),c- 1 (Si,...,s n+ i)- ( 37 ) 



Therefore 



Pr(S 1 ,...,S n+1 )=Pr(C- 1 (S 1 ,...,S n+1 )). (38) 



This implies that the probability distribution produced by taking the thermal ensemble for E'(S[, . . . , S' n+1 ) and 
applying C to the system is equivalent to a thermal ensemble with a new energy function 

E / (C- 1 (S 1 ,...,S n+1 )). (39) 



A quick calculation using Eq. (34) finds that 



A a " 

E'iC-^Sx, . . . , S n+1 )) = -(1 - St) + - £(1 - SiSi +1 ) (40) 

t=i 

which we see is equal to our original energy function E. Reversing the above argument we thus see that the energy 
function E can equally well be thought of as the ensemble corresponding to E' followed by the application of the 
mapping C . Further this later ensemble and computation has a clear interpretation in terms of a probabilistic 
initialization followed by probabilistic gates: when we apply controlled-nots with error target values then this causes 
a bit flip error in copying the information down the line. Thus we see that we can derive the same probabilistic 
gate expression for the one-dimensional Ising model as that which arose by reinterpreting the transfer matrix as a 
probabilistic circuit. 



III. THE PARTITION FUNCTION AS A PROBABILISTIC CIRCUIT 



We now turn to the more general setting of a circuit energy function for a generic circuit C = (G, L). We begin by 
focusing on circuit energy functions without a forced input. 

Let Ec be the energy function for the circuit C. For each logic gate, I, with inputs i±,i2,..-,ij and output 
ti, t-2, . . . , tk, define the following tensor, 

M (i) fl if Mh,i2,---,ij) = (ii,t 2 ,.--,*fc) / 41 n 

u,. | e -b otherwise ' ^ ' 

where b — j3A as before. Given this definition, we can now write the partition function Z as the contraction of a 
tensor network [31] plus a sum over all inputs and outputs. What is a tensor network? Take a graph and map a 
tensor to every vertex in such a manner that the tensor has an index for every edge of the relevant vertex. One 
can then perform a sum over two tensor indices connected via edges in the graph. A tensor network is thus a graph 
where vertices are tensors, and edges are tensor indices. Free edges are indices which have not been summed over and 
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connected edges are ones for which a sum has been performed. The entire network (graph) then represents itself a 
tensor with a number of unsummed indices equal to the number of free edges. 

Given the tensors M^> and the graph given by the circuit (which will have free edges for the inputs and outputs), 
we can build a tensor network out of these M^s. This tensor network will itself be a tensor network having indices 
for the inputs and outputs. Call this tensor i t t for inputs ii,...,i n and output t\,...,t m . Then it is 

easy to see that the partition function is equal to 

Z = E K,..,i n M,-,t m ( 42 ) 

ii,...,»„,ti,...,t m e{+i,-i} 

In other words, the partition function is given by the sum over all inputs, and all outputs of the value of the tensor 
network. 

What does this have to do with probabilistic circuits? Well now instead of using the tensors M, instead use 

p(l) _ n,...,lj,t u ...,t k , . 

ii,...,ij,ti,...,t k 1 + (2 fe — l)e~ fc 



or in other words 



J l + (2 fc ^l)e-'' if *2 j = (*l)*2j • ■ • ,tk) 



U < -' , ~-'» \ 1+{2 C 1)e - b otherwise 1 ' 

It is easy to check that this defines a stochastic matrix, i.e. a probabilistic gate. Indeed it is a probabilistic gate which 
performs the correct function of the gate with probability 1+ ^ 2 fcj_ 1 - )e -b and randomly flips the output to one of the 

-6 

other outputs with equal probability 1+ ^-i) e - b ■ 

Notice, importantly, that the difference between and M^> is a constant which depends on only on the number 
of output bits. So suppose we consider the tensor network for the circuit made now with P^s and not M^s. Then 
we can write the partition function as 

Zc = K J2 Piu...,i n ,tu...,t m (45) 

ii,;»„,ti,->tme{+i,-i} 

where K is the simple combinatorial factor 

K = + (2 k > - l)e- b ) (46) 
l 

and ki is the number of outputs of the Zth gate. Now examine Pi 1 ,...,i n ,t 1 ,...,t m - This is nothing more than the probability 
that we get output t\, . . . ,t m given that we had input i\, . . . ,i n to the probabilistic circuit with probabilistic gates 
P( l h In other words, we can express the partition function as 

Z C = J( ^ Pr(output|input) (47) 

inputs output 

Thus the partition function is really hiding a probabilistic computation. Notice further that we can explicitly calculate 
that partition function because output Pr(output | input) = 1, 

Z c = KT (48) 

where n is the number of input bits. 

Now we turn to the case where we force the input to the circuit. Let x\, . . . , x n be the inputs to the circuit C and 
E c .x be the energy function for the forced computation. Define the modified energy function 

Ec, x , y = E c , x {s)+E y {s), (49) 

where 



E y {s) = Y^nE Vi {s) (50) 
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with 

5SH- (51 » 

and the sum is over the output bits, ij. Call the partition function associated with this setup Zc, x ,y 

Note that Zc, x ,y\-y 1 =---=-y m =o is equal to the partition function for energy function with a forced input x, Zc <x - 
Define 

n 

Pin(ili ■ ■ ■ j in) = Pin,j(ij), (52) 



where 



if ij 



Then it is easy to see that 



PnnAij) = \ _£T if - J • (53) 

Z C ,x=-K" 51 Pi 1 ,...,i n ,t 1 ;-,t m Pin{il,---,in), (54) 

ti,;»„,ti,-.tme{+i,-i} 

where 

A"' = if (1 - e- b ) n . (55) 

Thus we see that we can interpret Zq, x as a sum over a probabilistic circuit, but now with a probabilistic input 

corresponding to each bit being flipped with a probability 1 ^ e _ b . Again, because we are summing over all outputs, 
we can explicitly perform this sum 

Z c , x = K'. (56) 

We are interested in computing the probability that the thermal ensemble for the forced input circuit has output 
given by t/i , . . . , y n . This can be calculated by 

1 d m Zc, X ,y 



Pr(h = yi,...,t m = y m ) = -- 

PZc,x,v 97i •• • din 



7i=-"=7m=0 



1 d m Z, 



C,x,y 



(57) 



71 = 



/3Z c ,x ■■ ■ dj n 
Note that 

Zc,x, v = K' e-^»Wp il ,..., iB , tl ,..., tm P in (ii,...,t n ) (58) 

»l,;*n,tl,-.tm6{ + l,-l} 

so that the partial derivatives may be evaluated only over the first term. There we find that 

Qm e -0E y (s) 



dji ■ ■ ■ d-/2 



'11^, (59) 



71 = 



Thus we find that 

Pr{h = yi, . . . ,t m = y m ) = ^ P Sl ,...,s n , yi ,...,y n Pin{si, ■ ■ ■ , s n ). (60) 

si ,. . . ,s n 

This is the main result of this section: the probability of the ground state spin computing model output bits is 
equal to the probability of running the probabilistic circuit described by P on inputs described by input probability 
distribution P in . The gates in P are simply noisy versions of the deterministic gates which fail with a fixed probability 
related to the temperature. 

This result is interesting in two manners. First it allows one to use techniques designed for probabilistically failing 
gates to be ported over to our model. For instance this will allow us to use methods for constructing fault-tolerant 
circuits. It is also interesting in that it shows that certain many-body quantum systems can be efficiently simulated. 
Computed circuit energy functions can be simulated by directly implementing the probabilistic gate that these systems 
correspond to. Here we have shown that output bits for the circuit energy function are related to the output bits 
for the probabilistic computation. A straightforward generalization shows that this is also true of other bits in the 
system. We can efficiently simulate systems with computed circuit energy functions by executing the probabilistic 
computation these energy functions represent. 
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IV. THE GATE MODEL DERIVATION 



Next let us turn to the same analysis as the prior subsection, but now using the gate trick as we did for the one- 
dimensional model in Section [II B| To do this we must first define the equivalent of the controlled- not gate. Consider 
a logical gate I with input i\, . . . , ij and output t\, . . . , tk which computes the function (t\, . . . , tk) = f(ii, • ■ ■ , ij)- 
Define the function C t from {0, l} j x {0, l} k to {0, lp x {0, l} fc as 

{{h, ■ ■ fi(h, ■ ■ if h = ••• = t k = 

(ii,...,ij-,0, ...,0) if (ti,...,tk) = f(h,...,ij) (61) 
(ii, . . . , ij,t%, . . . , tk) otherwise 

This function thus computes the function on input (ii , . . . , ij, 0, . . . , 0), uncomputes the function on input . . . , ij , fi(ii, 
and does nothing otherwise. Note that Ci is a bijection and is self-inverse, C, = C;. Defining an order to gates in 
our circuit such that the circuit computes properly, l\ first, h second, etc. Then we can define the function on all our 
bits of 

C c = C lr o...oC l2 o C h (62) 

Suppose we initially start with input bits initialized to an input i\ = x\, . . . , i n = x n and all other bits initialized to 
0. Then if we apply Cc to these bits we will place the result of the computation in their appropriate locations across 
the circuit. 

Next define an energy function on all of our system. This is most easily defined in terms of the inputs to the full 
circuit and the outputs to the gates / G L, 

E'{s')= A(^ 1+ ...+^ x J+E £ ^---'*m)- ( 63 ) 



and 

This energy function thus assigns for non-global inputs an energy penalty for output to gates being anything different 
that kl . For global inputs it adds a penalty for every input which is not the correct input from x. The thermal ensemble 
resulting from this system is then trivially described: the input qubits are in the correct input Xi with probability 
i+exp(-b) an d the internal output bits of a gate are grouped together and have a probability 1+ ( 2 k i-i) cxp(-b) being 

-b 

all zeros and probability zi) e xp(-b) *° ^ e an ything else. 

Imagine applying Cc to the state described by the ensemble for E'(s'). Clearly the ground state of E'(s') is such 
that the circuit C will be correctly computed. Further, because the ensemble related to E'(s') has erred inputs, the 
resulting probability distribution will have erred inputs. Finally the gates will function properly only when the output 
was properly initialized to all 0s. This occurs with probability 1+ ^ 2 k, CX p(_ b ) f° r logic gate I and otherwise the gate 

~ b 

fails by an equally probable error state with probability 1+ ^ 2 k, exp (_b) • I n other words the ensemble is exactly the 



one describing the probabilistic circuit in Section |III| Let us now show that this corresponds to the thermal ensemble 
of our original Eq construction. 



As in Section II B the trick here is to express the probability distribution for the ensemble after applying Cc as 

p r( s ) = Pr (s\s')Pr(s'), (65) 

s' 

where 

Pr(s\s') = S Cc (s'),s = s s ',c^(s)^ ( 66 ) 

such that 

Pr^^PriC^is)). (67) 

Thus the probability distribution resulting from after the application of Cq to the E'(s') thermal ensemble is equal 
to the probability distribution arising from E'iC^is)). 
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We will now show the E'(C is exactly the energy function we would define for a ground state spin computing 

construction for our circuit C . Recall that 



C'c 1 = C h oC h ---o C, 



It- 



(68) 



Note that in evaluating E'(C which is a sum of terms, we can evaluate the action action of Cq 1 on each of 

these terms separately and then re-sum. Thus we first examine 



Co 1 



E 



E 



AO^i,-.*! H h &i' m ,^x m ) 



(69) 



This is because the inputs to the circuit commute through the individual gate elements Ci i . Indeed this observa- 
tion along with the order of Cq 1 allows one to derive the term arising from each E[{Ci(ii, . . . ,ij t ,ti, . . . ,£&,)) term 
independently. In particular the only term which is not A for this energy function is 

E' l (C l (h,...,i jl ,f l (i 1 ,...,i jl ))) = (70) 

But this simply means that 

E' l {i 1 ,...,i jl ,h,...,t kl )=0 (71) 

only when (ti, .. = fi(ii, ■ ■ • , ij t ) and is A otherwise. This is just our normal definition of the circuit energy 

function. 

Thus we see that using the Cc construction one can construct a probabilistic circuit from the thermal ensemble for 
E' c x (s') which is the probabilistic circuit described in Section III Further we have shown that this can be thought of 
as arising from the original energy function Ec, x - This provides an alternative derivation of our main result. 



V. GENERALIZED CLASSICAL GROUND STATE SPIN COMPUTATION MODELS 



Having shown that classical ground state spin computing at non-zero temperature can be mapped to probabilistic 
circuits it is interesting to return to our original model and consider what generalizations of the energy function allow 
for this mapping to occur. The crucial assumption during our derivation of the probabilistic re- interpretation of the 
partition function was that we could convert the tensors of Eq. (41) to tensors of Eq. (43) in such a way 
that F" could be interpreted as a probablistic gate. This requires (1) the division of tensor indices into input and 
output indices and (2) for each input to the circuit the normalization to make the outputs probabilities (positive and 
sum to unity) is the same. 

Let us begin by showing that not all tensors and divisions into input and output vertices can result in a probabilistic 
interpretation. This can be demonstrated with a simple example on two bits: 



E( Sl ,s 2 ) 



if s% = s 2 
A if s% = s 2 
2A otherwise 



(72) 



If we chose Si as an input and s 2 as an output, this gives rise to the tensor 



M 



1 



„-26 



e -2b e -b 



(73) 



In order that the first column of this matrix sum to unity we must divide by 1 + e~ 2& . However when we do this the 
second column will not sum to 1. A similar argument occurs if we had chosen s± as the output and s 2 as the input. 
Thus there is no way to scale M in such a way that it can be interpreted as a probabilistic gate. 

A more relevant example of a failure for a model to not be amenable to our scaling arising in quantum-dot cellular 
automata models. In quantum dot cellular automata models bistable quantum dots are engineered in such a way that 
the ground state enacts a classical computation as in our model. Most of the treatments of these models deal with 
the quantum mechanics of these devices [9l [10l [32] • One can, however, build a classical model of these systems as was 
done, for example, by Wang in [llj . In this approximation one derives a classical statistical mechanical many spin 
system as in our model. When one does this for a quantum-dot cellular automata wire one ends up with exactly the 
model we can consider above in Section [llj However when one considers this for the majority gate constructions (see 
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[5HTT]) one obtains a model for which our technique cannot be applied. In particular the majority gate that they 
propose has a classical energy function which is 

E M (S 1 ,S 2) S 3 , S m ) = -A(5i5 m + S 2 S m + S 3 S m ) = -A(S 1 + S 2 + S 3 )S m . (74) 

Here Si,S 2 ,S 3 are the input spins and S m is the output spin. (Note that we have used one less spin than in the 
traditional presentation, but this does not affect our conclusions as one can think of the extra spin as being connected 
to the central spin by an identity gate.) Note that the lowest energy configuration for a given input correctly computes 
the majority function. That is if we fix Si, S2, and S3, then the lowest energy value of S m is given by the majority 
of the inputs, S rn — MAJ(S\, S2, S3). Note, however, that if we sum over all outputs for input Si = S2 = S3 = +1 

J2 exp(~bE M (S 1 ,S 2 ,S 3 ,S m )) = e- 3b + e- 3b , (75) 

S m |S , i=S 2 =S 3 =+l 

but if one uses input Si = S2 = +1 and S3 = — 1 one obtains 

exp(-bE M (S 1 ,S 2 , S 3 , S m )) = e b + e~ b . (76) 

S m \S 1 =S 2 =+l,S 3 =-l 

Thus we cannot properly globally normalize the entire tensor in such a way as to to obtain a probabilistic gate. 

However, the same two-spin energy function can be used to describe a fanout gate, in contrast to the many-spin 
interaction used in Eq. (JlJ, 

fe+i 

E F (S 1: S 2 , S k+1 ) = -ASi S i- ( 77 ) 

j=2 

Here Si is the input and the other spins are the output. By changing the direction of the circuit, we can show that 
this element can be represented as a probabilistic gate. I.e. the majority gate used in quantum dot cellular automata 
can not be mapped to a probabilistic circuit, but reversing the role of the inputs and outputs, this gate can be used 
to construct a fan-out gate. To see this note that, as for our previous fanout model, the ground state of this system, 
where the computation is performed, is degenerate. However now, as opposed to a single excited state there are 
multiple excited states. For example if we consider a fan-out to two spins we obtain an M as given in Tabled From 



Si 


s 2 


S3 


e -PE F 


+1 


+1 


+1 


e 2b 


+1 


+1 


-1 


1 


+1 


-1 


1 


1 


+1 


-1 


-1 


e~ 2b 


-1 


+1 


+1 


e- 2b 


-1 


+1 


-1 


1 


-1 


-1 


1 


1 


-1 


-1 


-1 


e 2b 



TABLE I. The 2 output fan-out Boltzman factor. 

this table we see that while the outputs that are in error have different Boltzman factors, we see that the sum over the 
output spins is indeed independent of the output. Indeed when we interpret this as a probabilistic gate by dividing 
out by K = 2 + e 2b + e~ 2b , then we see that this is a probabilistic gate which, instead of having a global failure where 
all possible output failures have equal probability, the fan-out spins have independent error probabilities. In other 
words it is as if the fan-out occurred and then with probability e - b+e b each of the outputs is independently flipped. 
We will use this fan-out gate in the next section when we discuss fault-tolerance. 

What is a necessary and sufficient condition for an energy function to allow for a proper normalization to a proba- 
bilistic gate? If we insist that this work for all temperatures, it is easy to give this answer. Let E^ . . . , ij, ti, . . . , £&) 
denote the energy function for a logic gate with inputs ix,...,ij and output ti, ...,tk- Let £(i±, . . . , ik) = 
{E(ii, . . . , ik,t%, ■ • ■ ,tk)\ti> ■ ■ ■ jtk S {0, 1}} denote the set of energies for a given input. Then the necessary and 
sufficient condition is that the £{i\, . . . , ifc)'s must be permutations of each other. The sufficiency of this condition is 
obvious. Necessity follows from noting that we are requiring this to hold for all temperatures and thus the sum over 
outputs of the exp(— j3E) must be a rearrangement of this sum and therefore a permutation of the output energies. 
Note that this condition implies that the ground states of a logic gate energy function must be degenerate for all 
possible inputs (but note that this alone is not sufficient for the condition to hold.) 
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A. Computational Complexity of Identifying Probabilistic Circuits in Thermal Systems 

Having established a necessary and sufficient condition for the possibility of scaling the M^ 1 ' tensors as probabilistic 
gates, a natural follow up question is how one can determine whether it is possible to take an energy function 
E(si, . . . , s n ) and determine whether or not it can be interpreted as a probabilistic circuit. Here we will show that 
this problem is computationally intractable, even when we are given information about the breakdown of E(s±, . . . , s n ) 
into terms which we wish to interpret as gates. 

Suppose that we are told the E(s\, . . . , s n ) is made up of a sum over logic gates and preparations 

E(s u ...,s n ) = J2 E u (78) 
i 

i.e. we are given a specification of the Ei in terms of a constant number of bits and we are told that each of these Ei will 
correspond to a logic gate or a preparation of initial inputs. We can then ask: can we interpret this as a probabilistic 
circuit with or without prepared inputs? Every Ei involves some inputs and output bits, but in general we will not 
be told which of these are inputs and which are outputs. In general an may thus have many different tuples 
of input bits and output bits such that we may interpret P™' as a probabilistic gate by proper global normalization. 
For example gates which have doubly stochastic P^' tensors may be run either backwards or forwards (a doubly 
stochastic matrix is a stochastic matrix whose rows and columns sum to 1.) Another example is given by a fan-out 
gate in our original energy function model 

E f ( Sl ,S 2 ,S 3 ) = {° A ^s 1 = s 2 = s 3 (7g) 

I A otherwise. 

In this case we may interpret the M^ L \Si, S 2 , S3) arising from this energy function as having input Si and outputs 
S 2 and 5*3, input S 2 and outputs S± and S3, or input S3 and outputs 52 and 53. Contrary to this, suppose that the 



energy function is that of the one-to-two fanout described in Eq. (77), 



E F (S 1 ,S 2 ,S 3 ) = -AS 1 {S 2 + S 3 ). (80) 

The M»(5i,S 2 ,5 3 ) for this energy function has only one possible orientation between inputs and output: Si is an 
input and S2 and £3 are outputs. This in general will lead to the following question: given the possible input and 
outputs tuples for the Aff's is it possible to assign these in such a way as to interpret the resulting total tensor as a 
probabilistic circuit. We will now show that this problem is NP-complete and is therefore (unless P = NP), in general, 
intractable. 

To see that the problem is NP-complete proceed as follows. Assign to each bit of our physical system a boolean 
variable, Xj. Assume that a given logical energy function E\ is involved in only a constant number of bits. Since Ei 
is involved in a constant number of bits we can explicitly calculate all of the sets of inputs and output bits for which 
can be globally normalized so as to make it a probabilistic gate. The variables x^s are going to represent whether 
the corresponding bit is associated is an input or and output. In particular, TRUE, will represent that the bit is an 
input and FALSE will represent that the bit is an output. Thus from the sets of input and output bits for which one 
can consistently label inputs to M« we can construct a boolean function on the relevant bits such that this function 
is true if and only if a valid assignment of boolean variables respects a possible input and outputs. Thus, for example, 



for Ef above in Eq. (79), the boolean function would be 

[(xi A ^x 2 A ^x 3 ) V (->x\ A ^x 2 A x 3 ) V (-1X1 A x 2 A -1X3)] (81) 

Constructing one such boolean function for each logic gate we can then construct a boolean function on all bits Xj 
which is the conjunction (logical and) of all of the gate boolean functions. Thus we see that by calculating for each E\ 
the possible combinations of inputs and outputs that have an allowed interpretation in terms of a probabilistic circuit, 
we have a one-to-one mapping of the problem to the problem of deciding whether there is a satisfying assignment 
for the boolean function made up of the conjunction of logical terms with only a constant number of involved bits. 
In order to demonstrate that the problem is NP-complete we need now only show that the boolean function we can 
construct in this manner are general enough so as to yield a version of the satisfaction problem which is NP-complete. 
In particular it is clear that this problem is in NP, since we can evaluate in polynomial time whether a given assignment 
of inputs and outputs can be interpreted as a probabilistic circuit. To demonstrate N P-completeness we need to show 
that it is at least as hard as another NP-complete problem. 

We will show that the problem is NP-complete by showing that some instance of the problem can be mapped to 
the one-in-three satisfiability problem of Schaefer |33| (now sometimes called the monotone one-in-three satisfiability 
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problem.) Define R(x,y,z) as the boolean function which is true iff exactly one of (x,y,z) is true. Notice that this 



is exactly the boolean formula of Eq. (81 1 (with xi = x, x^ = y, and x 3 — z.) In the one-in-three satisfiability 
problem one is given a boolean function which is the conjunction (logical and) of clauses each of which are of the form 
R(x, y, z) for a choice of variables x, y, z in the problem. Schaefer proved that the one-in-three satisfiability problem is 
NP-complete. Actually there is a slight subtlety here as the problem Schaefer considers must allow repeated variables. 
We will now show that it is possible to take a one-in-three satisfiability problem and convert it into an instance of the 
problem of spotting whether a given energy function can support interpretation as a classical probabilistic circuit. 

As mentioned above the one-in-three satisfiability problem has as input a conjunction of boolean functions which 
are all of the form R(x,y, z), and where some boolean variables can be repeated, 

r 

f(xi, . . . ,x n ) = f\ R(a,i,bi,Ci), (82) 



j 



where a.j, &j, Ci € {xi, . . . , x n } and none of the dj, b iy and c,; are the same variable for fixed i. Let = {aj 
Xi} U {bj\bj = Xi} U {cj|cj = Xi} denote the set of a, 6, c variables that are Xi. We will define an energy function on 
3r + X)i=i 1-^*1 spins which, when converted int o the problem of trying to assign consistent inputs and outputs will 



yield exactly the boolean function / of Eq. (82). Call the first 3r spins Si and the last X to t = Yli=i 1^1 s P ms ^i- 
First define the energy function 

r 

Em(Si, ■ ■ ■ , S 3r ) = ^2 Ef(S3j-2> Ssj-ii S 3 j), (83) 
3=1 



where Ef is the fan-out energy function defined in Eq. (82). When one converts this into a logical expression for 
whether the spins are inputs or outputs, one gets exactly the one-in-three logical clause of Eq. pll). 



However, at this point, the variables that will arise in these clauses are all different. This is what the other X tot 
spins are for. Define the following energy function for the "identity" gate: 

E?\p 1 ,P 2 ,...,P k ,Q 1 ,Q 2 ,...,Qk) = {\ ifPl=Ql,P2 = Q2 '""" ,Pfc = Qfc • (84) 

I A otherwise 

This energy function corresponds to a logical reversible gate which acts as identity from the P, spins to the Qi spins. 
This energy function has the property that the Pj variables are all either inputs or all outputs (and the Qi variables 
are also all either inputs or outputs, opposite to the Pj variable assignment.) Thus we can use gates of this form to 
effectively make the spins of the Em gates constructed above the same variable. In particular define the following 
energy function 



E e (Si,..., S 3r ,Ti,. . . ,T Xtot ) — '^2, X * (ftsQti^iJXip], • • • , S x,[\x t \}-,T x (i) ,T x (i) +v . . . ,T X ^ + ^ X ^), (85) 



where — J2j<i 1^1 an d Xi[j] is the jth element of Xi translated over to a Si variable, i.e. flj corresponds to 
S 3 i-2, bi corresponds to S 3 i-\ and Cj corresponds to S 3 i. 
Now consider the energy function 

E(Si, . . . , S 3r , Ti, . . . , Tx tot ) = E e (Si, . . . , S 3r , Ti, . . . , Tx tot )) + Em (Si ... , S 3r ), (86) 

and stipulate that every term in the sum of E e will correspond to one gate and every term in Em will correspond to 
one gate. Then if we convert this problem over into a boolean satisfaction problem as describe above, where each spin 
is a boolean variable representing whether the spin is an input or an output, we will obtain a boolean expression that is 



equivalent to the problem described in Eq. ( 82 ) . Thus we have shown that we can, for a given one-in-three satisfaction 
problem, Eq. (82) construct an energy function, along with labeling of the terms in this function corresponding to 
different gates, such that determining whether this energy function can be reinterpreted as a probabilistic computation 
is equivalent to the original one-in-three satisfaction problem. This implies that the problem of determining whether 
such probabilistic computations exist is at least as hard as this NP-complete problem. Combined with the fact the 
problem is in NP this means that the problem of deciding whether a given energy function can be thought of as 
enacting a probabilistic circuit is N P-complete. 

Thus we have shown that the problem of determining a many-spin statistical mechanics system enacts a ground 
state spin computation, even given the decomposition of this system into logical gates, is NP-complete. We note in 
passing that this result is of perhaps some philosophical significance: determining whether a system can be thought 
of as a computer is shown to be computationally intractable (assuming P 7^ NP). 
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B. Models that have a Probabilistic Gate Interpretation 



So far we have considered spin models in a very general setting unrelated to physical theories. We point out, here, 
however, that there are a variety of models that have arisen in physics which admit such an interpretation. We have 
already seen that that one dimensional Ising model admits an interpretation as a series of identity gates which err 
with some probability. Here we point out some other models which admit such interpretations. 



1. Ising Models on Trees 

A tree is a connected graph, T = (V, E), with vertex set V, edge set E, and no cycles. Assign to each spin a vertex, 
S„ for v £ V. To each edge e = (v, w) assign an non-zero energy for the Ising coupling, J : E — > M. — {0}. Then the 
Ising model on the tree has an energy function given by 

E(si, . . . , s\v\) = ^ J( v ,w)S v S w . (87) 

Fix a root vertex v r € V. First consider the model at zero temperature. At T = 0, following our discussion of the 
fan-out in Section [Vj it is easy to see that the ground state of this model can be thought of as a series of fan-out 
gates originating from the root vertex v r with possible bit flips applied to the fan-out depending on the sign of J( v . w y 
At T 7^ we see that, as in the discussion of the fan-out gate, each of the individual fan-out wires (with possible bit 
flips), will fail independently with a probability of 

= e pj M+e -pj M - (88) 

Note that a failure for an antiferromagnetic coupling (bit flip) corresponds to an identity gate. Ising models on 
trees have been studied in a variety of settings [Ml |3S] and in fact the interpretation of this model in terms of a 
probabilistic circuit has been used extensively in the computer science literature, where the model goes under the 
name of broadcasting |12j . 



2. Zi Lattice Pure Gauge Theories 

Perhaps slightly more interesting that Ising models on trees is the fact that lattice gauge theories [351 [37] can 
be cast as a probabilistic gate circuits. For definiteness, first consider a (pure) Z% lattice gauge theory on a two 
dimensional square lattice with closed boundaries. In this model one places spins on the edges of the lattice and 
defines the energy function 

E=J2 J P SiS 2 S 3 S4 (89) 

pep Si,s 2 ,s 3 ,Si.eSp 

where P is the set of all plaquettes of the square lattice, and Sp are the spins surrounding plaquette p. Consider 
an individual term in this sum, and assume for simplicity that for all plaquettes, p € P that J p — J < 0. We can 
view the individual term JS1S2S3S4 as a gate element with either 0, 1, 2, or 3 inputs and 4, 3, 2, and 1 outputs 
respectively. For example if we view it as a gate with 2 inputs and 2 outputs, then, in its ground state this corresponds 
to a probabilistic gate which, with equal probability, changes the input to a set of spins with the same parity of —Is 
as the input. In other words the logical evolution is 

(+1, +1) -> 50% (+1, +1) and 50% (-1, -1) 
(+1, -1) -> 50% (+1, -1) and 50% (-1, +1) 
(-1, +1) -> 50% (-1, +1) and 50% (-1, +1) 

(-1, -1) -> 50% (+1, +1) and 50% (-1, -1). (90) 

The gate with 3 inputs and 1 outputs is deterministic and is given by /(Si, S2, S3) = S1S2S3. The gate with inputs 
corresponds to a preparation of 4 spins which have an even number of —1 spins and this state is prepared with equal 
probability. The gate with 1 input and 3 outputs is a probabilistic gate which outputs an even number of —Is if 
the input is +1 all such possibilities occurring with equal probability. Similarly it outputs an odd number of —Is if 
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the input is —1 all such possibilities occurring with equal probability. By picking boundary terms to serve as input 
and output bits, it is clear that it is possible to order the energy terms in the above sum in such a way that we can 
interpret the ground state as a probabilistic cellular automata like model using the above gates. At finite gate these 
probabilistic gates all fail with some probability, where failure results in an equally likely failed output for the gate. 
While we have define this for the two dimensional square lattice, it is clear that the above construction can result 
in a probabilistic circuit for a far greater number of lattices, the main criteria being that there exists a way to label 
inputs and outputs in such a way that a proper circuit is constructed. 



VI. MAKING CLASSICAL GROUND STATE SPIN COMPUTATION ROBUST AT FINITE 

TEMPERATURE 



Having shown how the computed circuit energy functions give rise to thermal ensembles which can be interpreted 
as probabilistic circuits with probabilistic inputs, we now discuss how it is possible to make the model of classical 
ground state spin computation fault-tolerant at finite temperature. Indeed this is now rather simple given that we 
have a mapping between the behavior of these systems at finite temperature and probabilistic circuits which have 
gates and preparations failing with a fixed probability. We will need, however, to use the more general class of energy 
functions considered in Section IVl 

The general theory of computing in the presence of gates which fail goes back at least as far as the work of von 
Neumann [8]. Von Neumann considered a model in which each logic gate failed with exactly the probability e. In 
order to get this model to compute reliably, von Neumann used two major techniques. The first techniques is that 
one needs to suitably encode information: in order to do this one encodes a single bit across a bundle of wires in a 
logic circuit and interprets this as (or 1) when a majority of the wires are in the (or 1.) Given such an encoding it 
is then easy to compute in a parallel manner such that the effect of gates failing is simply to decrease the proportion 
of Os (or Is) in an encoded (or 1.) However, the fact that the gates fail does cause the ratio of correct bits in an 
encoded bundle of bits to decrease. In order to deal with this von Neumann used a second technique whereby faulty 
gates were used to perform error correction. Actually this portion of von Neumann's analysis is not quite satisfactory, 
as it requires the use of a random permutation. However Pippenger |38j has shown that one can de-randomize this 
construction (one can obtain reasonable parameters for this method by using the expanders defined in 39 ) . The end 
result of this is that one can show that if one wishes to compute a logical circuit of size n to an accuracy 1 — 5 (that 
is, have the circuit fail with probability 8) with gates that fail at e < e% for some fixed threshold e-t, one can do this 
using a circuit with the unreliable gates which is only a 0(log fc 4) larger, where A: is a fixed constant. This effectively 
means that one can create for all effective purposes reliable logic gates out unreliable logic gates with an overhead 
which is logarithmic in the inverse of the error desired. 

Having described von Neumann's construction one can now see how it is possible to make classical ground state 
spin computing fault-tolerant. In particular we have shown above how the thermal ensemble for such systems gives 
rise to probabilistic gates. For gates constructed using energy functions like Eq. ([I]) we have shown that the thermal 

(2 k — l)e~ b 

ensemble arising for these functions is equivalent to a circuit in which gates fail with probability 1 _j_( 2 fc -i)e-'> ^ e 

gate has k output bits. For k = 1 this is ^j-^=h ■ We would like to use this for von Neumann's e failure probability. One 
complication arises, however, which is that in von Neumann's model the fan-out gates do not fail. These gates are 
needed during the error correcting stage of his procedure. We get around this by using fan-out gates constructed from 
energy functions like those in Section [V] In particular define the fan-out gate via Eq. 
this gives rise in our probabilistic gate setting in which each fan-out resulting wire fails with probability e -b +e b ■ 
These independent failures can then be reinterpreted as errors for the gates that follow them in von Neumann's error 
correcting construction. One should note that there is a subtlety here in that if we desire to recover von Neuman's 
error model where the gates fail with exactly probability e then a complicated calculation is needed in order to figure 
out e in this model. However, for large enough j3 the effect will mostly be that the error probabilities add such that e 
will be the sum of the fan-out error and the gate error. 

Putting this together we can conclude that there is a critical temperature below which one can use fault-tolerant 
constructions to design ground state spin energy functions which compute correctly with probability 1 — S by using 
0(log fe 1) more interactions in the energy functions. This means that while the general model of ground state spin 
computing is not tolerant to errors, for a cleverly designed setup one can overcome this difficulty and build robust 
devices, assuming that the temperature is low enough to keep the error probability below the error threshold. 



77 1 . As shown in Section M 
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VII. RELATIONSHIP TO COMPUTATIONAL COMPLEXITY THEORY 



Now we turn to a relationship between our model and computational complexity. This connection is provided 
because the ground state spin computing model at zero temperature is intimately related to the Cook-Levin theorem 
from computational complexity. 



The Cook-Levin theorem [22J 122] is a fundamental result in the theory of computational complexity. Here we briefly 
review this result. For more details one is referred, for example, to the classic textbook |40) . 




Traditionally the theory of computational complexity is founded upon the computational model of a Turing ma- 
chine [41] . The Church- Turing thesis roughly states that the Turing machine model adequately captures the notion of 
what can be performed by any physical instantiation of a computer and decades of experience have yet to challenge 
this assumption. A deterministic Turing machine consists of the machine and a linear tape onto which elements from 
a finite alphabet, T, can be written into different cells on the tape. The machine itself has an finite set of internal 
states, Q, a method for reading a symbol written in a cell on the tape, a method to write a symbol in a cell on the 
tape, and the ability to move either left or right one cell. Formally the tape is assumed to be infinite or at least 
infinitely extendable. Initially an input, consisting of a string of symbols from a finite alphabet £ C T is placed on 
the tape and the physical machine is placed on the first symbol of this string. The Turing machine is assumed to start 
itself with its internal states in a special state called the start state, qo € Q. A deterministic Turing machine then 
proceeds to "compute" as follows. At each time step the deterministic Turing machine examines the symbol of the 
cell it currently occupies and its internal state. Then conditional on these two variables the Turing machine writes a 
symbol onto cell (possibly from a larger, but still finite, alphabet than the alphabet of the initial strings) it occupies, 
moves either left or right, and modifies its internal state. These rules can be specified by a function, the transition 
function, 6:QxT^QxTx {L, R}, where L and R denote moving left or right. The deterministic Turing machine 
does this indefinitely until it reaches one of a set of its internal states known as halt states. These halt states are 
labeled as cither "reject" or "accept." If wc initially input the string w and the deterministic Turing machine halts 
in accept state, then we say the deterministic Turing machine accepts w. Similarly if on input w the machine halts 
in a reject state, then we say that the machine rejects w. Note that on some strings the Turing machine may neither 
accept or reject a string, but continue indefinitely. The language of a Turing machine is the set of all strings taken 
from the input alphabet for which the Turing machine accepts for that string being input. 

Deterministic Turing machines capture many important features of real modern computers. For example, the 
running time of most algorithms on a variety of hardware are closely related to the similar problems cast in terms of 
deterministic Turing machines. An important class of languages are those for which a deterministic Turing machine 
accepts a string from the language in a time bounded by a constant times a polynomial in the size of the string. Such 
languages are said to be in the complexity class P. 

In contrast to a deterministic Turing machine, non-deterministic Turing machines are a model of computing which 
is not thought to properly encapsulate what can be efficiently computed on a modern computer. A non-determinisitic 
Turing machine works in a manner similar to a deterministic Turing machine: it has an internal state, it can read 
and write symbols onto a tape, etc. The difference is that at each step in the computation, a non-deterministic 
Turing machine doesn't have to chose one new state, new symbol and direction to move, but may branch to multiple 
such triples. The transition function is now not S : Q x V — > Q x T x {L, R} but instead is of the form 6 : Q x T — > 
P(Q x r x {L, R}) where P denotes the power set. Non-deterministic Turing machines thus explore multiple pathways 
at once. Such a Turing machine halts when one of the states of any of the multiple branches it is exploring enters 
a halt state. Again we can define whether the initial input is accepted or rejected depending on whether the state 
where the computation halts is a reject or accept state. Given this we can define the language of the non-determinstic 
Turing machine as the set of strings in which the machine accepts the given input. The time complexity of a non- 
deterministic Turing machine is the number of steps of executed by the Turing machine irrespective of the number of 
branches explored by the machine. 

Because a non-deterministic Turing machine is allowed to explore multiple pathways in its computation it is not clear 
that such machines can be reasonably mapped (in terms of say, time complexity of problems) to modern computers. 
This question can be put on a more formal footing by defining the complexity class NP. NP is the class of languages 
which can be recognized by a non-deterministic Turing machine in a time bounded by a constant times a polynomial 
in the size of the string being decided. The statement of whether NP is strictly larger than P is the famous P = NP 
problem, one of the most important open problems in the theory of computational complexity. That being said, most 
researchers do believe that P 7^ N P and that the model of a non-deterministic Turing machine is more powerful than 
computers we can physically construct. 



A. 



The Cook-Levin Theorem 
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Variable 


Represents 


H ik 

Sqk 


True if tape cell — p(n) < i < p(n) contains symbol j at time < k < p(n) 

True if the Turing Machine is at tape cell — p(n) < i < p(n) at time < k < p(n) 

True if Turing machine is in state q at time < k < p(n) 



TABLE II. Variables in the Cook-Levin theorem. 



Given the above brief overview of the basic models of deterministic and non-deterministic Turing machines we can 
now introduce the Cook-Levin theorem. To define the Cook-Levin problem we introduce the SAT problem. 

Boolean Satisfaction Problem (SAT): Given as input a description of a Boolean function from n 
boolean inputs to one output, / : {0, 1}™ — > {0, 1}, decide whether there exists a set of inputs x\, X2, ■ ■ ■ , x n , 
with each Xi e {0, 1}, such that f(xi,x 2 , ■ ■ ■ , x n ) = 1. 

Cast in terms of languages, the SAT language is the set of strings describing a Boolean function such that the Boolean 
function has some input which evaluates to 1. The Cook-Levin theorem then states that the SAT is complete for NP. 
This means that every language in NP can be converted by a deterministic Turing machine in polynomial time to a 
SAT language and every SAT language can be converted into a language in NP. Thus, in a real sense, SAT captures 
what it means for a lanuage to be in NP. For example, if we could efficiently solve a SAT problem in polynomial 
time we would be able to solve every problem in NP in polynomial time and vice versa. This is especially important 
considering that numerous problems beside SAT have been identified as NP-complctc. This class of problems are 
essentially "the most difficult" in NP and are correspondingly the problems for which no known polynomial time 
algorithms for all instances of the problems in the class are thought to exist. 

We are most interested in the proof of the Cook-Levin theorem. Before reviewing this it is useful to provide another 
equivalent characterization of NP which uses, instead of a non-deterministic Turing machine, a prover and a verifier. 
The prover is assumed to have unlimited computational power while the verifier is limited to running a polynomial 
time algorithm on a deterministic Turing machine. A language, L, is in NP if the prover can supply to a verifier a 
polynomial sized proof p that a instance x is in the language. This proof must always be correct, i.e. the verifier must 
always be convinced, by using a Turing machine from P, that the instance x is the in the language. Similarly the 
prover must never be able to convince the verifier incorrectly. In other words if x is not in L, then the verifier must 
always reject. For example SAT is in NP by this definition since the prover can simply supply a satisfying input to the 
boolean function if the function can be satisfied, and the verifier can compute the boolean function on this input in 
polynomial time. Similarly if the function does not admit a satisfying assignment the verifier can never be convinced 
that the instance is not in the language because no string the prover can give will lead to a satisfied boolean function. 
Note that in the prover-verifier setting we do not require that the prover try to convince the verifier that an instance 
x is not in L, only that if it is true that x is not in L that the prover cannot convince the verifier that x is in L. 

Given the prover/verifier definition of NP we can now discuss the Cook-Levin theorem. The Cook-Levin theorem 
states that SAT is NP-complete. A problem, P, is NP-complete if and only if every problem in NP can be transformed 
into P by a deterministic Turing machine in polynomial time. Thus SAT, since it is NP-complete, captures much of 
what is difficult about problems in all of NP: if there were a polynomial time algorithm for SAT then there would be 
a polynomial time algorithm for every problem in NP. 

The proof that SAT is NP-complete has two directions. We have already stated the first direction, that SAT is 
a NP problem. The meat of the proof is in the other direction, that every problem in NP can be converted into 
a SAT problem. Here we do this in a slightly non-standard manner, using the prover-verifier definition of NP. In 
particular for any problem, L, in NP we know that there exists a verifier who takes as input the problem instance, x 
and a certificate p, runs these on a deterministic Turing machine and verifies whether the instance is in the language. 
Imagine now the time history of such a verification by a Turing machine. Initially the Turing machine is at the 
beginning of the string, in its initial state, and the input string is on the tape. At the next time step, the machine 
can change its state, move left, right or stay in place, and possibly change the symbol on the tape where the machine 
is located. Notice, however that the contents of a cell on the tape, whether the machine is at that location, and the 
internal state of the machine if it is on that cell at the next time step is a function only of the contents of the cell 
and its neighbors at the prior step, whether the machine was on one of those cells, and its internal state if it was. 
In other words the computation can be described locally in space and time. Given this we can now describe how to 
convert a verifier, which is just a deterministic Turing machine and the strings representing the instance x and the 
proof certificate p into a problem in SAT. 

Suppose that the verifiers deterministic Turing machine runs in time p(n) where n = \x\ + \p\ is the size of the 
instance and the proof and p is a polynomial function. Define the following boolean variables, j fe, Hi.k? Sqk to 
represent the history of the deterministic Turing machine in verifying the proof. We think about the first two variables 
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as being arranged in a large 0(p(n)) x (p(n) + 1) tableau onto which the computational history of the deterministic 
Turing machine is carried out. We can now define a SAT problem over these variables which is satisfiable iff the 
deterministic Turing machine accepts the input string x,p. This SAT formula is conjunction (logical and) of a clauses 
which enforce (a) that the initial contents of the tape are x, (b) that the machine is in the correct initial state, (c) the 
initial location of the head of the tape, (d) that there is only one symbol per tape cell, one state per time, one tape 
position per time, (e) that the tape cell is unchanged unless it is the one which can be written upon, (f) the machine 
properly updates according to the Turing machine transition function, and (g) the machine finishes in an accepting 
state. For example condition (b) consists of clauses with a single variable S q0t o where qo is the initial state of the 
Turing machine. The most complicated clause in this construction is the one that encodes the local update rule for 
the machine, (f). This can be done by adding the clauses 



for valid transitions of the Turing machine (that is for q,q' , j, f ,d representing a valid change of state, symbol, and 
direction of movement for the read/write head of the Turing machine.) Putting this together one thus sees that once 
can take the deterministic verifier along with the instance and x and construct a polynomial sized boolean formula / 
which is satisfied iff the verifier accepts the proof on the instance. Note that we have not forced a constraint on the 
proof input to the circuit, but that the definition of NP in a provcr- verifier setting leads to the satisfiability of the 
formula / iff the instance x is in the language. 

We can now see that connection between the Cook-Levin theorem and the classical ground state spin computing 
model. In particular a central part of the proof of this theorem is that one can construct a satisfiability formula 
directly from the history of a computation. In a similar manner classical ground state spin computing constructs 
an energy function whose ground state is the history of a computation. The conversion of a logic circuit to such an 
energy function is essentially the arithmetization technique of computational complexity [421 143) . Note however that 
the important role here of the input and outputs to the constructed formula: in the case where the output is forced, 
we are led to NP-complete problems, but when the input is forced then we are led to efficient constructions. Motivated 
by this it is interesting to consider how the above constructions works but now at finite temperature. 



We will now show that it is possible to use classical ground state spin computing to define, in analogy with 3 — SAT 
a Promise-MA-complete problem. MA stands for Merlin-Arthur and is essentially a probabilistic generalization of 
the complexity class NP [201 HI]- It is defined as the follows: 

Promise-MA: A language L = L yes U L no C £*, L yes n L no — 0, is a promise problem in Merlin- Arthur 
(Promise- MA) if there is a probabilistic polynomial-time Turing machine V and a polynomial p such that 
for all strings x € L, 

• If x € L yes then there exists a w such that \w\ — poly(\x\) and V(x, w) accepts with probability at 
least 2/3. 

• If x 6 L no then for all w such that |iu| — poly{\x\) then V(x, w) accepts with probability at most 1/3. 

Here E is the alphabet over which the language is defined. This version of the definition of MA is a promise problem, 
since L yes U L no does not necessarily cover all possible elements of E*. The reason that this class is called Merlin- 
Arthur is because it can be defined in terms of an interactive proof system between a computationally unbounded 
prover, Merlin, and a more limited verifier, Arthur, who can perform feasible probabilistic computations (polynomial 
time probabilistic computations with bounded error, BPP.) Merlin is trying to convince Arthur that a certain instance 
a: is in a language L yes . He does this by communicating a polynomial size proof w to Arthur. Arthur can then take 
this proof and run a polynomial time probabilistic computation on this proof the result of which is: if a; € L yes , then 
Arthur is, with high probability, convinced the x is in L yes , while if x e L no no matter what proof Merlin supplies 
Arthur cannot be convinced that x € L yes except with a small probability. 
We define the following problem: 

Boundary Expectation Value (BEV): Let k be a positive integer and f3 a k digit number. Suppose 
we are given an energy function on n bits, E(s\, . . . , s n ) which can be written as a sum of energy functions 



(i?i,fc A S q ^k A iij.fc) — > (ifj+cz,fc+i A S q '^k+i A Lij> : k+i) 



(91) 



B. A Promise-MA Problem Prom Classical Ground State Spin Computing 




(92) 
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in such a way that this energy function can be interpreted as a circuit energy function (see Eq. Q) with 
j > input bits and 1 output bit, and each term in Ei is specified with k bits of precision. Consider the 
thermal ensemble generated by this energy function at inverse temperature (3. We promise that either 
there exists input bits i%, ia, . . . , ij and output bit t such that for the reduced probability distribution on 



these bits 



or for any input bits ■ ■ ■ ,ij, 



Pr(ii,i 2 ,...,i J -,t = l)>| (93) 



Pr(i 1 ,i 2) ...,i J) t = l)<l (94) 



The problem BEV is to determine whether the first of these conditions, Eq. (93), holds. In other words 
yes instances of this problem satisfy Eq. (93) for some inputs and no instances satisfy Eq. (94) 

for all i%, . . . , ij. 

We claim that BEV is Promise-MA-complete. The proof of this is fairly straightforward given our prior discussion 
of the Cook-Levin theorem. 

First we claim that BEV is in Promise- MA. To show this we must show that there exists an interactive proof of 
the form described above for the problem BEV. To see this note that a computationally unbounded Merlin can sample 
from the thermal ensemble for an instance of this problem and can check whether there exists an input ij, . . . , ij such 
that Eq. (93) holds or whether for all inputs i\, . . . ,ij, Eq. (94) holds. This will then serve as the proof, w, in the 



definition of MA. Given the description of (3 and the circuit energy function, Arthur can construct a probabilistic 
Turing machine which implements the circuit related to this energy function using our construction of a probabilistic 
circuit from Section [III| Then Arthur can run this probabilistic Turing machine on the proof that Merlin supplies. It 
is clear that if a; is a yes instance of BEV, then the proof w Merlin supplies will convince Arthur that he has a yes 
instance of this interactive protocol. Further if a; is a no instance of BEV it is also clear that no proof that Merlin 
supplies will convince Arthur that X is a yes instance. Thus it is clear that BEV is in Promise-MA. 

The other direction of the proof requires that we show that every problem in Promise-MA can be converted into 
a BEV problem. To do this we follow the idea of the Cook-Levin theorem in that, for a MAproblem x £ L yes , we 
convert the probabilistic Turing machine V that verifies the proof w of the MAproblem into a thermal circuit for 
which the associated BEV problem has (x,w) as the satisfying input. Similarly, for an MAproblem x € L no , the 
thermal circuit corresponding to V must correspond to a no instance of BEV, since otherwise there would exist 
inputs w = ii,i2, . . . ,ij which would cause V(x,w) to accept with probability greater than 1/3. Clearly to do this 
we will make use of our probabilistic circuit representation result from Section |III| The main tool that we need in 
such a construction is that we need to be able to construct circuit energy functions which perform probabilistic gates 
of a form we desire, for a fixed inverse temperature /3. To do this we can proceed as follows. We can use our fault- 
tolerant constructions to produce circuits which operate with near deterministic behavior. Then in order to augment 
a deterministic computation we need a source of randomness, indeed in most standard definitions of a probabilistic 
Turing machine one requires randomness which is close to uniformly random across the random bits being used (or 
in other words one requires bits which are nearly equal probability of being in or 1.) To obtain such randomness in 
our constructions, notice that if we do not force our inputs to be or 1 then there is, in a thermal ensemble, an equal 
probability to be and 1. This allows one to create random bits in an circuit energy function. However in order to 
get this to work with the deterministic part of our Turing machine, we must find a method to get this randomness 
into encoded bits in our circuits. This can be accomplished by simply taking the equal probability and 1 bits and 
running them through the circuit which performs error correction on the bundle corresponding to these bits. Thus 
we see that we can construct at finite /3 and circuit energy function which is near deterministic and which can also 
have input bits which are nearly uniformly random. Thus we can proceed just as in the Cook-Levin theorem. For 
a tableau describing the history of the probabilistic Turing machine, V, (when this probablistic Turing machine is 
simply a deterministic Turing machine aided by bits of randomness, we can construct an energy circuit function which 
implements the probabilistic Turing machine spatially across this tableau. 

Thus we have shown that BEV is Promise-MA-complete. While the proof of this result was rather straightforward 
and followed quickly from understanding the Cook-Levin theorem it is interesting to note that very few MA-complete 
problems are known. In fact the only other non-natural problem which is Promise-MA-complete to our knowledge 
is the stoquastic 6-SATproblem The nearest comparable result is the classic result for the complexity of finding 
the ground state and computing the partition function for the Ising model of Barahona |45j . 
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VIII. QUANTUM MODELS 

Finally we would like to briefly mention the quantum models which show some similarity with our classical model. 
In particular these models were the inspiration for investigating this problem and thus this serves as a good location 
to introduce the open problems concerning these problem. 

One of the quantum models relevant to this discussion are universal adiabatic quantum computing schemes. In these 
models one shows how to adiabatically drag a many-body quantum system from one easily preparable ground state 
to the ground state of another many-body quantum system whose ground state is a superposition over the history of 
a quantum circuit. Thus, similar to our models, the end result is a ground state which encodes a computation: but in 
this case the computation is a superposition over the history of the computation. These models show that adiabatic 
quantum computing, which previously had been only used for optimization problems |46H48| can also be used for 
universal quantum computation [T5HI5] , However the final state of these models are systems with small energy gaps 
and thus the effects of working at non-zero temperature will have a profound effect on these models. Indeed it is for 
exactly this reason that such universal adiabatic quantum computing schemes are not known to be fault-tolerant (see, 
however, |49j.) Similar reasoning follows for an earlier model of ground state quantum computing due to Mizel and 
co-workers [T31 IT4"] . 

The problem of having a system whose ground state correctly computes but whose excited states do not is exactly the 
problem we have addressed in this paper for the classical ground state spin computing. An interesting question then, 
when considering the quantum models, is whether there is a similar reinterpretation quantum thermal ensembles 
as spatially enacted quantum computations. This is an important open problem which might conceivably lead to 
fault-tolerant methods for adiabatic quantum computing. 

IX. CONCLUSION 

We have introduce a set of energy functions on many bits (spins) which have the property that their ground state 
can be thought of as a spatially distributed deterministic computation. At temperature greater than zero we have 
shown that the thermal ensemble arising from these models can be reinterpreted as a spatially distributed probabilistic 
computation. Further, above zero temperature we see that the gates of a ground state spin computer can become 
unreliable and fail to execute the desired computation with high-fidelity. However with the mapping of the thermal 
ensemble to probabilistic circuits we have shown how it is possible to make versions of the desired deterministic circuits 
which are fault-tolerant. We have shown that determining whether a given classical energy function can be thought 
of as enacting a ground state computation is NP-complete. Finally we have shown that a problem concerning the 
thermal ensembles arising in our model give rise to a rare complete problem for the complexity class Promise- MA. 
The models we have considered here are mostly devoid of connections to actual physical systems of interest, besides 
connections to quantum-dot cellular automata. An interesting an important open question about these models is 
whether they arise in naturally occurring physical systems, or a suitably engineered system. Another interesting 
question is the rate at which the ground state spin model thermalizes: a proof that the system reaches thermal 
equilibrium in a time polynomial in the size of the circuit being implement would be another step towards making 
this model more physically relevant. 
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